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A GENERAL METHOD FOR DETERMINING COINCIDENCE 
CORRECTION OF COUNTING INSTRUMENTS 


By T. P. Kohman 


ABSTRACT 

The method herein described for determining and applying coincidence corrections in counters 
is an extension of the general method of paired sources. Several pairs of radioactive sources are 
used, and for each pair, counts are taken with the instrument on the sources separately and together. 
A general relationship between the recorded counting rate R and the “true” rate N is assumed, 
being the first few terms of an infinite power series: 

N = R+ r R 2 + yR 3 + <pR 4 + .... 

where ris the resolving time of the counter and v, <p, . . . . are independent parameters. By a 
least-sqiares method, the values of the parameters for the instrument in question are determined 
which most adequately correlate the experimental data. Subsequent measurements made with the 
instrument can then be corrected by means of this equation. 

A few examples of the use of this method are presented. They indicate that the method is a 
satisfactory solution of the problem and should be of considerable value. 

INTRODUCTION 

All instruments which count random events, such as those encountered in the measurement of 
radioactivity, have finite resolving powers, and hence are unable to separately distinguish and record 
events occurring very close together in time. As a consequence some counts are missed, the fraction 
lost increasing as the counting rate increases. This results in a nonlinear variation of the response 
of the instrument to the intensity of the source of events. The difference between the recorded rate 
and the “true” rate (that is, the rate which would be observed if the instrument could count all events 
to which it is normally sensitive) is called the coincidence loss or correction. 

In practice coincidence losses are often appreciable, and therefore, in order to obtain accurate 
measurements in radioactivity with counting instruments, it is necessary to correct the recorded 
counting rates for this effect. Moreover, the useful upper limit of a counter is greatly extended by 
a knowledge of its coincidence corrections. Coincidence loss measurements also supply valuable 
information about the mode of operation of counters, and provide a useful guide in improving in¬ 
struments so as to decrease their losses. 

The problem has received attention for a number of years, and several methods of determining 
and applying coincidence corrections are described in the literature. However, all of these methods 
are either inconvenient, not generally applicable, or not sufficiently accurate. Recent experience at 
the Hanford Engineer Works and other Metallurgical Project Sites has indicated the importance of 
the problem and the inadequacy of the previously described methods. This report describes a new 
method which is capable of any desired accuracy and which is almost completely general in its 
applicability. 
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THEORY OF THE METHOD 
Notation 

To facilitate the discussion, the notation to be used is presented herewith: 

R = recorded rate. 

N = true or corrected rate (including background), 
b = background rate. 

N' = N - b = corrected rate due to source alone, 
t = time of count. 

Subscripts A and B refer to individual members of a pair of sources. 

Subscript C refers to the combined sources of the pair. 

r = resolving time of counter. 

v,<p , . . . . = supplementary parameters. 

Methods Involving the Use of Paired Sources 

Reviews of methods of determining coincidence corrections have been given by Beers 1 and 
Kohman. 2 The use of paired sources, first proposed by Moon, 3 has been frequently employed. 4 * 5 * 1 * 2 * 8 * 7 
If each of two radioactive sources of approximately equal strength is measured separately, and 
then the two are measured together, taking care to keep geometrical factors constant throughout, 
the combined counting rate (corrected for background) will be less than the sum of the individual 
rates (corrected for background), since the fractional counting loss is greater at the higher 
counting rate. If a simple mathematical relationship between recorded and true rates involving 
but a single parameter applies for the instrument, a single set of measurements on one pair of 
sources can be used to evaluate the parameter. 4 * 1 * 2 * 6 * 8 * 7 However, it has been shown 2 * 8 that the 
simple one-parameter relationships previously proposed do not adequately express the behavior 
of real counters. 

In order to avoid this difficulty, a completely general method which involves no assumed 
mathematical relationship between true and recorded rates was developed. 2 This method em¬ 
ployes a number of paired sources of different intensities covering the range of the instrument, 
and the data are analyzed by a graphical method of successive approximations which gives a 
curve of fractional correction versus recorded rate. Unfortunately, this method is tedious, time 
consuming, and possibly subject to personal factors, A completely mathematical analysis of the 
data would be a considerable improvement. A prerequisite of this is a mathematical relationship 
between true and recorded rates. 

Single-parameter relationships between N and R which have been proposed are these: 

a) Johnson and Street; 9 Ruark and Brammer: 10 

N = R + r NR (1) 

b) Volz, 11 Schiff: 12 

N = Re TN ® 

c) An approximation of (b) due to Hull: 4 

N = R + tR 2 + 1/2t 2 R s 


(3) 
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d) An approximation of (a) 6, V: 

N = R + tR 2 


(4) 


For each of these relationships, methods have been developed for evaluating the parameter t 
from a set of paired source measurements. The derivations of equations 5 to 9 proceed from the 
fact that N' A + N' B =N'o 


a) Beers: 1 


R A + R b - R C - b [(R a + R B - R C - b)(Rc - b)/(R A - b)(R B - b)] 2 
2 (R a - b) (R B - b) 8(R C - b) 


Kohman: 8 



Equation 5 is 


->P"rJ^ (R A + RB-Rc-b)J (6) 

more exact at low counting rates, and equation 6 at high rates. 


b) Crawford: 6 

N c = (R A+ R B ) 2 /Rc (7) 


t is then obtained from N^, R c , and equation 2. For equation 7 to be 
applicable, the background must be small enough to be neglected, and R A 
and Rjg must be very nearly equal. 

c) Hull: 4 


1/2 (R a s + Rg 3 -R c 3 )r 2 + (R A 2 + R b 2 - R C 2 ) T + R A + R B — R c — b = 0 (8) 


t is obtained by solution of this equation by the quadratic formula, 

d) Jarrett; 8 Kohman; 2 * 3 Metcalf and Hennessee: 7 

r = (Ra + R b - R c - b )/< R c 2 - «a 2 " r b 2 > 0) 


Inspection of equations 1, 2, 3, and 4 reveals that they are all nearly equivalent for low counting 
rates, where the difference between N and R is very small (~1%). Moreover, any one will adequately 
represent the behavior of actual counters in the region of small losses. Equation 4 means that the 
fractional correction is proportional to R, and because of the simplicity of this equation and its 
corollary 9, these equations have received extensive use in rough evaluations and applications of 
coincidence corrections. 6 * 13 * 8 ’ 7 In general, however, no single-parameter equation, theoretical or 
empirical, can adequately relate N and R when the losses are large. As an illustration, a semi- 
theoretical equation which expressed the observed behavior of a Geiger-Mueller counter over a 
wide range of counting rates contained three parameters (equation 2). However, that equation is 
not suitable for practical application, nor is it applicable to other types of instruments. 

In theory, the relationship between two dependent variables such as N and R can be expressed 
by an infinite power series: 


N = c 0 + c x R + c 2 R 2 + C 3 R 3 + c 4 R 4 + • • • • 


(10) 
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If we set c 0 = 0 
Ci =1 
c 2 = r 

equation 10 becomes equivalent to equations 1, 2, 3, and 4 for small values of R. Hence, we are 
justified in rewriting the series as follows: 

N = R + tR 2 +uR s +4>R 4 + • • • . (11) 

For practical purposes, the number of terms of equation 11 to be retained will be determined 
by the behavior of the instrument and precision required. 

As previously mentioned, when only two terms of equation 11 are used, giving equation 4, the 
parameter can be evaluated from measurements on a single pair of sources by means of equation 9. 

It has been pointed out by Jarrett 8 that one might similarly use the first three terms of equation 11 
and evaluate the two parameters from measurements on two pairs of sources. One obtains two 
equations analogous to equation 9 but involving the two unknown parameters, and their simultaneous 
solution yields the values of the latter. Similarly, this method can be extended to any number of 
terms of equation 11 by employing as many pairs of sources as there are parameters. 

The Least-squares Evaluation of Paired Source Data 

The use of equations such as equation 9, or its higher analogues according to the method of 
Jarrett, has a serious disadvantage in that only a limited number of measurements can be admitted 
into the evaluation of the parameters. It is generally desirable to employ a greater number of 
measurements in obtaining the final result, so that an error in one measurement will not seriously 
effect the values of the parameters (as would be the case when the minimum number of measure¬ 
ments is used), and so that the entire counting range of the instrument can be represented in the 
data. The method of least squares provides a means of doing this. 

Consider a pair of sources which have been counted separately and together. Let us define the 
quantity 6, which we will call the residual, thus: 

6 =N' a + N' b -N' c 

If the N values were the true values, as would be the case if the correction of the R values 
were perfect, then 6 would equal zero (except for statistical fluctuations due to the random dis¬ 
tribution of counts). Actually, when the N values are obtained by an equation such as will be used 
in practice, the corrections will not be perfect, and 6 ^ 0. 8 is the algebraic sum of the errors in 
the corrections to R A , Rb> ^ R C* K ma y be considered as a measure of the error in the correction 
to R c since this is the largest correction, and hence, is expected to have the largest error. Simi¬ 
larly, we may define the fractional residual £ =6 /Rq, which may be considered a measure of the 

fractional error in the correction to R^. 

If we have a number of such paired source measurements, the best values of the parameters 
in equation 11 are those which make X £ 3 a minimum: 

a 

dr 

a 

dv 

sV 

1 


2 £a=0 

2 £2=0 
^ £ 2 = 0, etc. 
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One-parameter Case 

For the case where only two terms of equation 11 are retained, giving equation 4, the evaluation 
of the parameter t by the least-squares method proceeds as follows: 

£ j N'a + N'b - N'c 
R C ” Rc 

_ N a + N b -N £ . r _b 
R C 

RA+RR-Rr-b-t- t(Ra 2 + Rr 2 -Rr 2 -b) 

Rc 


D + rE 
Rc 


where D = Ra + Rb ™ Rc ~ b 
E =R A 2 + Rb 2 -Rc 2 


= (D 2 + 2tDE + t 2 E 2 ) 




E 2 


R C 

= G + 2tH + t 2 L 


where 


G = 


H 


2Lt R c 

~Zj r c 2 

_ y» E 2 

Rp : 


The condition for best fit is 


dr 


'y e 2 = 0 = 2H + 2 tL 


= -(h/l) 


whence: 


( 12 ) 
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Two-parameter Case 

For the case where three terms of equation 11 are used: 

N = R + tR 3 + v R s 

the least-squares evaluation of the two parameters follows: 


(13) 


r a + r b _ Rc _ b + t ( r a j + r b* _ R c’> + v (Ra* + r b’ - R c s) 
R c 

D + tE + V F 


*c 

where D = + Rg "Rc 

E = R A * + R B , -R C * 

F =R A * + R B , -H C S 

2 , ’2s?® 4,E * ,,n ‘ 


+ 2tDE + 2uDF + 2 t» EF +r 2 E a +u 2 F 2 ) 


m 


= G + 2tH + 2u J + 2 tu K + t 2 !, + u 2 M 


where 


‘£$ 

" 2 & 
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The condition for best fit is: 


d Ve 2 = 0 = 2H + 2uK + 2 tL 
dr -SLi 

_s_ y € 2 = 0 = 2J + 2TK+ 2oM 
dv 

Solving these simultaneous equation yields: 


r = (JK - HM)/(LM - K 2 ) 


u = (HK - JL )/(LM - K 2 ) 


(14) 


Three-parameter Case 

The derivation of the method for the case where three parameters are used 


N =R +TR 2 + uR 3 +<j>R 4 (15) 

proceeds similarly. One obtains three simultaneous equations in three unknowns, whose solution 
yields the desired values. The solution is straightforward, but the computation is somewhat longer 
than in the two-parameter case, especially since the solution of simultaneous equations containing 
more than two unknowns is somewhat involved. Consequently, the two-parameter case is probably 
the most satisfactory for routine practical use. Only if one needs high accuracy and the experimental 
measurements have been made with a compatible degree of precision, or if the fractional corrections 
are quite large (^50%) should it be necessary to invoke additional parameters. Obviously the method 
can be extended to any number of parameters. 

Significance of the Parameters 

It is customary to consider the coincidence losses of a counter in terms of quantities called the 
resolving time, the insensitive time, and the dead time. The resolving time is defined as the minimum 
interval between two events which can both be registered by the instrument. The insensitive time is 
the interval following a recorded event during which the instrument is incapable of recording another 
event. The dead time is a quantity peculiar to Geiger-Mueller counters, and is the interval following 
a discharge during which the tube is incapable of producing an electron avalanche. 

The insensitive time is not a property of the instrument, but of a particular event recorded by 
the instrument. With practically all types of counters, this quantity is not a constant, but may be 
different for different events. For example, the insensitive time of a Geiger-Mueller counter dis¬ 
charge may be shortened by the occurrence of another discharge a short time previously, 14 > 2 > 6 or 
it may be lengthened by the occurrence of an unrecorded discharge near the end of the interval. The 
insensitive time of a pulse in an alpha particle counter may depend on the length and position of the 
track of the particle. If the insensitive time had a constant value, the resolving time would have this 
same value. But, in view of the variability of insensitive times, the simple definition for resolving 
time given has no clear meaning. Nevertheless, it is desirable to retain the idea of the resolving 
time as a characteristic of the instrument, in the sense that the insensitive time is a characteristic 
of an individual event. 

Equations 1 and 2 are derived theoretically from certain physical assumptions, among which is 
that of a constant insensitive time or resolving time,r. The parameter has the same meaning in other 
approximate equations, such as equations 3 and 4. Since the empirical equation 11 becomes essen¬ 
tially equivalent to these equations at very low counting rates, we may call t in the latter also the 
resolving time. Physically, it is the limiting average value of the insensitive time as the counting 
rate approaches zero. With Geiger-Mueller counters, the insensitive time is constant at very low 
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rates, but with proportional counters or counting ionization chambers, the insensitive time varies, 
so the average is specified. Mathematically, the resolving time may be defined as 


d(N/R) 

dR 


R = 0 = — 


d(R/N) 

. dN . 


N =0 


The parameters v of equation 11 have no simple physical significance, but are merely 

constants employed to correlate the experimental observations. 

The parameter r, being the resolving time, has the dimension time and is always positive. The 
parameters v ,<£, .... have the dimensions time 2 , time 3 , . . . and may have either sign. 

The discussion thus far has been confined to counting instruments. The method is equally ap¬ 
plicable to some other types of measuring instruments which do not count. All radiation measuring 
instruments are expected to deviate increasingly from linearly as the radiation intensity increases. 
Examples to which the method can be applied are Lauritsen quartz-fiber electroscopes, high pres¬ 
sure gamma-ray ionization chambers, beta-ray ionization chamber, DC amplifier assemblies, 
photocells, etc. In such cases, the notion of a resolving time no longer applies, and it is more 
logical to rewrite equation 11 as 


I = M + aM 2 + bM 3 + cM 4 +. (16) 

where M and I are the measured and true intensities respectively. 

Once the parameters have been evaluated for a particular instrument by the method described, 
subsequent measurements made with the instrument can be corrected by the use of the corresponding 
equations 4, 13, 15.The coincidence correction should be added before subtracting the back¬ 

ground, especially if both are large. To facilitate the application of corrections, one may employ a 
graph or table showing the correction as a function of the recorded rate. 


APPLICATION OF THE METHOD 
Procedure for Taking Data 

The experimental part of the determination of coincidence corrections of an instrument consists 
of taking a series of sets of paired source measurements. Each set is comprised of four measure¬ 
ments involving two approximately equal sources: R^j Rb* Rc k. The series should contain three 
or more sets, with counting rates covering the range of the instrument for which the corrections are 
desired. 

The procedure for taking a set of measurements is discussed. Sources and counting conditions 
are chosen which will give the desired counting rates. With many instruments the counting rate can 
be adjusted by varying the distance between source and detector, or by inserting suitable absorbers. 

In order to eliminate possible changes due to variations in geometry, a dummy source should occupy 
the place of each active source when the latter is not in position. The order of measurements is: 

1) Place both dummies in position. Measure b. 

2) Replace dummy A by active source A. Measure R^* 

3) Replace dummy B by source B, with care not to disturb the position of source A. 

Measure R q\ 

4) Replace source A by dummy A, with care not to disturb the position of source B. 

Measure Rg. 
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The sequence is illustrated in the figure below, where S A and S 3 are the active sources and 
Da and Dg the dummies: 



I: b H : R a III: R c IV : R B 


It is generally desirable to count R A , Rb» and Rc * or e( l ua ^ periods of time. In this manner, one 
can obtain ihe greatest statistical accuracy for a given total number of counts . 1 This is desirable 
with instruments such as Geiger-Mueller counters whose lifetime is a function of usage. The great¬ 
est statistical accuracy for a given total expenditure of time is obtained when: 

*A : : k: : tb = ^5 a : ■ \f R C ■ ^« 1 *• 1: 

However, the loss in accuracy by making t A = tg = tc is slight, so even for instruments which 
are not damaged by usage, the convenience of using equal counting periods is sufficient to justify it 
It is also desirable to make tb greater than indicated by this ratio, since a more exact value of th. 
background is a useful indication of whether the instrument is functioning properly. 

The counting time necessary for each set depends on the accuracy required in the corrections, 
and may be estimated as discussed. Consider an instrument for which equation 4 holds, and 7 has a 
value of approximately 10 ~ 5 minute. We wish to determine r from a single set of paired source 
measurements using equation 9, the counting rates being approximately 1000, 1000, and 2000 per 
minute. We require an accuracy (fractional probable error) of 0.5% in corrected measurements at 
2000 c/m. The coincidence correction at this rate is approximately 40 c/m, and this quantity must 
be known with an accuracy of 10 c/m, or 25% of its value. This means that 7 must be known with an 
accuracy of 2, which in turn means that D ( = R A + Rr “ r C ~ b) must also have an accuracy of 
25%. Since D~20 c/m, its probable error must be 5 c/m. The probable error in D is equal to 


0.6745 




5a. 

tA 


+ 5el + 5c. 

t B tc 


tb 


^ 0.6745 


\[¥ 


when the counting times are equal. From this it follows that t = 73 minutes. The general formula 
for the counting interval is: 


8 K*__ _ 3.64_ 
P 2 R C ' P 2 ~R C 


(17) 


where t = t A - tg = tc = required counting interval 

P - fractional probable error required at R - Rc 

K = 0.6745 = factor used in computing probable errors. 

(If it is desired to express P in terms of some other measure of 
deviation than the probable error, the corresponding value of K 
must be used instead.) 

It is to be noted that the time requirement is independent of the resolving time. From equation 17, 
Table 1 has been prepared. 

In cases where equation 4 does not hold, and the method described herein is to be used instead, we 
may imagine each set of paired source measurements as establishing a point on the correction curve. 
Hence, the same considerations apply, and equation 17 and Table 1 can be used as a general guide in 
paired source measurements. Thus, when a series of such measurements is made, each set should 
contain approximately the same total number of counts. 


10] 
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Table 1. Length of counting interval t in minutes required for a given 
fractional probable error P for various combined source counting rates Re 
in county/minute. 


Rc\p 

5% 

2 % 

i% 

0.5% 

0 .2% 

0 .1% 

1,000 

1.5 

9.1 

36 

146 

910 

3640 

2,000 

0.7 

4.6 

18 

73 

455 

1820 

5,000 

0.3 

1.8 

7.3 

29 

182 

728 

10,000 

0.15 

0.9 

3.6 

15 

91 

364 

20,000 

0.07 

0.5 

1.8 

7.3 

46 

182 

50,000 

0.03 

0.2 

0.7 

2.9 

18 

73 

100,000 

0.015 

0.09 

0.4 

1.5 

9.1 

36 

200,000 

0.007 

0.05 

0.2 

0.7 

4.6 

18 

500,000 

0.003 

0.02 

0.07 

0.3 

1.8 

7.3 

1 ,000,000 

0.001 

0.01 

0.04 

0.15 

0.9 

3.6 


Table 2. Paired source measurements on Geiger-Mueller counter. Argon alcohol- 
filled, bell type, mica window tube with self-quenching circuit and Offner scaler. 


No. 

R a (c/m) 

Rg (c/m) 

R c (c/m) 

b 

t jl (min) 

5 (c/m) 

€ (%) 

I 

2334 

2485 

4750 

30 

3.6 x 10- 6 

- 6 

-0.13 

n 

3077 

3155 

6133 

28 

3.9 x 10- 6 

- 4 

-0.07 

in 

3089 

3116 

6134 

28 

2.3 x 10- 6 

—33 

-0.54 

IV 

3727 

3696 

7293 

33 

3.8 x 10- 6 

- 9 

-0.12 

V 

3779 

3692 

7324 

29 

4.6 x 10-» 

+ 12 

+ 0.16 

VI 

5452 

5356 

10505 

33 

5.2 x 10- 6 

+ 56 

+ 0.53 

vn 

5334 

5411 

10508 

34 

3.9 x 10 

-14 

-0.13 


Result of computation: r = 4.12 x 10 ” 6 min, and N = R + tR 2 . 
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Table 3. Paired source measurements on alpha counter. Air-filled, parallel-plate 
ionization chamber with linear pulse amplifier and Offner scaler. 


No. 

R a (c/m) 

R b (c/m) 

Rc (c/m) 

b 

Tj (min) 

6 (c/m) 

«(%) 

I 

294.4 

313.4 

603.4 

2.0 

13 x 10- 6 

+0.6 

+0.10 

II 

546.5 

538.2 

1076.6 

2.0 

11 x 10- 6 

+0.4 

+0.04 

IH 

736 

725 

1450 

2 

9 x 10 * 6 

-2 

-0.14 

IV 

1323 

1143 

2436 

2 

9.7 x 10- 6 

-18 

-0.74 

V 

2459 

2514 

4782 

2 

17.9 x 10 - 6 

+13 

+0.27 

VI 

4454 

4327 

8131 

2 

23.5 x 10- 6 

+1 

+0.01 

VII 

5924 

6245 

10860 

2 

29.8 x lO' 6 

-7 

-0.06 

Result of computation: r = 8.20 x 10“ 8 min, u 

= 1.147 x IO" 9 min 2 , and N = R + 

tR 2 + uR ! 

Table 4. Paired source measurements on alpha counter. Nitrogen-filled, parallel-plate 
ionization chamber with high-frequency amplifier and Offner scaler.* 

No. 

Ra (c/m) 

R b (c/m) 

RC (c/m) 

b 

t i (min) 

6 (c/m) 

«(%) 

I 

41,287 

41,111 

81,094 

2 

4.10 x 10' 7 

+70 

+ 0.09 

II 

80,308 

70,495 

145,492 

3 

5.44 x 10- 7 

-135 

-0.09 

HI 

108,794 

123,160 

216,621 

2 

7.73 x 10- 7 

+74 

+ 0.03 


Result of computation (see Table 6): r = 1.814 x 10" 7 min, v = 1.662 x 10 ~ 1Z min 2 , and 
N = R + tR 2 + v R 3 


♦ The losses in this instrument occur chiefly in the scaler, not in the chamber or 
amplifier. 

Table 5. Paired source measurements on Geiger-Mueller counter. 4 Hydrogen-filled, 
thin-walled glass tube (Eck and Krebs), with Neher-Harper quenching circuit and 


scaler 

. Department of Chemistry, 

University of Wisconsin, June 1942. 


No. 

RA (c/m) 

Rb (c/m) 

R c (c/m) 

b 

7i (min) 

6 (c/m) 

*(%) 

I 

579 

647 

201 

16 

13 

x 

10- 6 

+1 

+0.06 

n 

2,054 

2,491 

4,386 

18 

11.0 

X 

io - 6 

+39 

+0.88 

HI 

5,560 

7,090 

11,960 

20 

10.8 

X 

lO" 6 

+12 

+0.10 

IV 

7,320 

9,970 

16,120 

20 

10.6 

X 

io - 6 

+66 

+0.41 

V 

17,880 

18,460 

33,100 

20 

7.4 

X 

io- 6 

-470 

-1.42 

VI 

28,800 

32,600 

53,700 

- 

6.3 

X 

10- 6 

+630 

+1.18 

VII 

53,500 

58,200 

93,200 

- 

7.6 

X 

10- 6 

-430 

-0.47 

vni 

65,700 

77,400 

114,000 

- 

10.8 

X 

10 

+200 

+0.18 


Result of computation: t = 1.217 x 10” 5 min, u = 8.76 x 10“ n min 2 , <f> - 5.29 x 10 18 min 3 , 

and N = R + TR 2 + n R 3 + (j>R 4 . 
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The choice of counting rates for the various sets of paired source measurements can be made 
from the following considerations. The R c values should be fairly evenly distributed over the useful 
range of the instrument. The maximum useful rate of a counter for general purposes may be con¬ 
sidered to be that for which the coincidence loss is about 10 or 20 per cent. Some Geiger-Mueller 
counters, however, have lower useful limits because high rates may permanently alter the tube 
characteristics. Since sets with low values of Rc require long counting intervals and enter into 
the computations with less weight than the higher sets, it is usually not profitable to use Rc values 
lower than that counting rate for which the coincidence correction is about 1 or 2%. If more than 
3 or 4 sets are taken, duplication is permissible and has the advantage of indicating the consistency 
and reproducibility of the data. 

Sample Data 

To illustrate the application of the method to actual counting instruments, several series of 
data taken with various types of instruments are presented in Tables 2, 3, 4, and 5 together with 
the results of the computations. The examples have been chosen to illustrate the one-, two-, and 
three-parameter cases. 

Preliminary Examination of Data 

In order to obtain an idea of the consistency of the data and to determine the number of param¬ 
eters which should be computed, a preliminary examination of the data should be made by applying 
equation 9 to each set of measurements. The value of r so obtained is designated as n> since it 
implies a one-parameter relationship, equation 4. The values of t x are included in Tables 2 to 5. 

If, as in Table 2, the variations in r 2 show no trend with counting rate, and are no greater than 
expected from statistical fluctuations, a one-parameter computation should suffice. The data of 
Tables 3 and 4 indicate a progressive increase in t 1 (the initial decrease in Table 3 is not signifi¬ 
cant), indicating that two parameters should be computed. The data of Table 5 show a significant 
decrease in followed by a significant increase, suggesting that three parameters are probably 
required to express the behavior of the counter. 

Sample Computation 

To demonstrate the manner in which the computations are set up and carried out, a typical 
two-parameter computation is shown in Table 6. The computations involving a different number of 
parameters follow similar patterns, except that with three or more parameters the solution of the 
simultaneous equations must be handled differently. 

In the computation, sufficient significant figures should be carried to allow for the loss which 
occurs in solving the simultaneous equations. In the three-parameter computation of Table 5, 
three significant figures were lost in the solution of the equations, so eight figures were carried 
up to this point in order to give satisfactory accuracy in the result. Because of its length, this 
computation was consequently quite tedious. 

Checking the Computations 

The computations should be checked against the data by computing the residual (6) and the 
fractional residual (e =5 /Hq) for each set. If the values of e are all small^that is, of the order 
of magnitude of P in equation 17 and are distributed randomly in sign, it is assumed that the 
equation obtained satisfactorily represents the behavior of the instrument. If, however, the 
values of € are all of the same sign, or show a significant trend, an error in the computations 
is indicated, or an additional parameter must be computed. 

In Tables 2, 3, 4, and 5, there are given for each set of measurements the values of 5 and £ 
computed from the result of the least-squares computation. In each case, these quantities are small 
and of random sign, indicating that the equation obtained is adequate for the correction of measure¬ 
ments made with the instrument. 
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Table 6. Sample computation, illustrating method of laying out work for deter¬ 
mining two parameters with aid of calculating machine. Data from Table 4. 


Set No. 

I 

n 

in 

r a 

41,287 

80,308 

108,794 

r B 

41,111 

70,495 

123,160 

R c 

81,094 

145,492 

216,621 

b 

2 

3 

2 

R A + Rb ” r C “k = D 

1,302 

5,308 

15,331 

r a 2 

1.70461 x 10 9 

6.44937 x 10 9 

1.18361 x 10 10 

r b 2 

1.69011 x 10 9 

4.96954 x 10 9 

1.51684 x 10 10 

R c 2 

6.57623 x 10 9 

2.11679 x 10 10 

4.69246 x 10 10 

r a 2 + r b 2 ~ r c 2 = E 

-3.18151 x 10 9 

-9.74901 x 10 9 

-1.99201 x 10 10 

r a 3 

7.03782 x 10 13 

5.17936 x 10 14 

1.28770 x 10 15 

r b s 

6.94821 x 10 19 

3.50328 x 10 14 

1.86814 x 10 15 

R c 3 

5.33293 x 10 14 

3.07976 x 10 15 

1.01649 x 10 16 

r a 3 + r b s - r c 3 = f 

-3.93432 x 10 14 

-2.21150 x 10 15 

-7.00903 x 10 15 

de/r c 2 = V 

-6.29960 x 10 2 

-2.44463 x 10 s 

-6.50821 x 10 s 

DF/R c a = W 

-7.78970 x 10 7 

-5.54549 x 10 s 

-2.28996 x 10 9 

ef/r c 2 = X 

1.90338 x 10 14 

1.01852 x 10 15 

2.97542 x 10 15 

E 2 /R C 2 = Y 

1.53918 x 10 9 

4.48996 x 19® 

8.45636 x 10 9 

f 2 /r c 2 = z 

2.35376 x 10 19 

2.31044 x 10 20 

1.04692 x 10 21 

Calculation of data in Table 6: 



X V = H 

-9.58280 x 10 s 

HK 

-4.00971 x 10 19 

£ W= J 

-2.92240 x 10 9 

JL 

-4.23324 x 10 19 

£ X = K 

4.18428 x 10 15 

HK-JL 

+2.2353 x 10 18 

a 

it 

>* 

1.44855 x 10 10 

LM 

+1.88529 x 10 31 

£ Z - M 

1.30150 x 10 21 

K 2 

+1.75082 x 10 31 



LM-K 2 

+1.3447 x 10 30 

JK 

- 1.2281 x 10 25 

JK-HM 

+1.814 x 10' 7 

HM 

-1.24720 x 10 25 

LM-K 2 7 


JK-HM 

+2.439 x 10 2S 

HK-JL 

LM-K 2 “ V 

+1.662 x 10' 12 
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Figure 1. Coincidence corrections for alpha counter. Data of Table 3. 

A: Found by present method. B: Previously assumed for this type of instrument. 

Graphical Illustration of Results 

For two of the examples considered, the fractional correction has been plotted as a function of 
the recorded rate. 

Figure 1 shows the behavior of the alpha counter from which the data of Table 3 were taken. For 
comparison, the dash line indicates the previously assumed losses of such instruments, 0.8% per 
1000 counts/minute. 13 

The difference between the two curves indicates the necessity of determining coincidence cor¬ 
rections separately for every instrument. 

In Figure 2 the continuous line shows the results of the three-parameter computation based on 
the Geiger-Mueller counter data of Table 5. The dash line is the result of the independent analysis 
of the same data by the method of successive graphical approximations. 2 The two curves coincide 
within the accuracy of the measurements, indicating the validity of both methods. The dash curve 
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RECORDED RATE (c/m) 

Figure 2. Coincidence corrections for GM counter. Data of Table 5. 

A: Found by present method. B: Found by method of successive graphical approximations. 


is perhaps the better of the two since it yields a value for L € 2 of 0.00005, whereas that obtained by 
the present method is 0.00046. A fourth parameter would presumably reduce the difference 
between the two curves. 

To facilitate the correction of counter data, it is often convenient to use a graph showing the 
correction (N-R) as a function of the recorded rate R, instead of applying the correction equation 
to the individual results. The plot may be made on Cartesian or logarithmic coordinate paper, 
depending on the range of values. Such plots are shown in Figures 3 and 4. 

AC KNOWLEDGMENTS 

The data of Tables 2 and 3 were obtained through the cooperation of Robert Jolly. The computa¬ 
tions were made by Robert R. Jones, Miss Virginia Webb, and Mrs. Gertrude Paas, through the 
cooperation of Benjamin Butler. 



CORRECTION (c/m) 


— i— a 


BIIKIBII 

H 


00 2 


Figure 3, 





CORRECTION (c/m) 


MDDC - 905 


l 



RECORDED RATE (c/m) 


Figure 4. Coincidence correction chart for GM counter, logarithmic plot. Data of Table 5, 
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